function deg = rad2deg(rad)
deg = rad/pi*180;